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Abstract. We present minor axis kinematic profiles for a well-studied sample of 19 early- to intermediate-type disk 
galaxies. We introduce, for the first time, the use of single-burst stellar population (SSP) models to obtain stellar 
velocities, velocity dispersions and higher order Gauss-Hermite moments (/13, ^4) from galaxy spectra i n the near- 
infrared Can triplet region. SSP models, which employs the synthetic spectra of Vazdek fs et al.l {2003), provide a 
means to address the template-mismatch problem, and are shown to provide as good or better fits as traditional 
stellar templates. We anticipate the technique to be o f particular use for h igh-redshift galaxy kinematics. We 
give the measurement of a recently defined CaT* index (|Cenarro et al.ll2001af) . and describe the global properties 
of the bulge kinematics as derived from the kinematic profiles. We detect small-amplitude minor-axis rotation, 
generally due to inner isophotal twists as a result of slightly triaxial bulges or misaligned inner disks; such inner 
features do not show peculiar colors or distinct CaT* index values. Velocity dispersion profiles, which extend well 
into the disk region, show a wide range of slopes. Flattened bulges tend to have shallower velocity dispersion 
profiles. The inferred similarity of bulge and disk radial velocity dispersions supports the interpretation of these 
bulges as thickened disks. 

Key words. Galaxies: bulges - Galaxies: kinematics and dynamics - Galaxies: structure - Methods: data analysis 



1. Introduction 

The measurement of kinematics in galaxies has been one 
of the major achievements in astrophysics in the second 
half of the last century. Our knowledge in this field has 
evolved from cru de kinematic determinations from pho- 
tographic plates (jBurbidge fc Burbidgd ll96Sh. to high- 
resolution, detailed measurements of the li ne-of-sight ve- 
locity distributions (hereafter LOSVDs; i.e. iBender et al.l 
Il99i lFisheJll997t IVega Beltran et al.ll200lh . 

In the seventies, methods were devised to measure 
velocity dispersio ns in addition to ra dial velocities, i.e. 
cross-correlation ([Tonrv fc Davisl ir979). Fourier Quotient 
([Sargent et al.lll977l) . In the past decade, with the im- 
proved quality of the data, a new set of algorithms 
were developed allowi ng us to meas ure LOSVDs (Fourier 
Correlation Quotient, Bender 1990), and quantify devia- 
tions from pure single Gaussian ([van der Marel fc Franxl 
Il99fl and multiple Gaussian LOSVD w ith UGD 
([Kuiiken fc Merrifieldl 11993 IPrada et al.l ll996). Galactic 
kinematics, however, suffers from an important technical 



constraint: template mismatching. Stars are observed to- 
gether with the object which kinematics we are interested 
in, to be used as spectral templates. However, it is dif- 
ficult to choose the right star, since galaxy light comes 
from a wide range of stellar types. Although, for ellipti- 
cals in the V-band, it has been found that h alf of the light 
comes from G and K giants (Pickles 1985) and, there- 
fore, those stellar types are good templates, this situa- 
tion may change dramatically from galaxy to galaxy or 
even within a galaxy. Therefore kinematic parameters will 
be sensitive to the stellar template used. Much work has 
been devoted to address the templat e mismatch problem 
(iRix fc White! Il992t iGonzale A 119931 IStatler et all Il996t 
IStatler fc Smecker-Handll999l) . 

In the last decade, ne w evolutionary stellar pop- 
ulation synthesis models (IBruzual A. fc Chariot) 1993; 
[Worthev et"all|l994[ iBuzzonil Il995t iTantalo et all Il996t 
IVazdekis et al.lll996[) have helped us to significantly in- 
crease our galaxy stellar populations knowledge. These 
models predict line strengths of a limited number of 
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spe ctral features, such as those of the Lick/I DS system 
fi.e lWorthev et al.l Il99i IVazdekis et all Il99rll) . However 
since these models were unable to predict full SEDs 
(Spectral Energy Distributions) at high enough resolu- 
tion, they could not be used as templates for galaxy stellar 
kinematics. This situation has recently changed with the 
development of new SSP model SE Ds in the blue an d vi- 
sual at resolution 1.8 A (FWHM) bv lVazdekisl dl999h and 
in the near-IR Ca H triplet regio n at 1.5 A (ly azdekis^LaL 
2003). In fact. IVazdekisl {l999) and IVazdekis fc Arimoto 
(1999) have shown the potential use of the SSP models as 
templates for determining the main kinematic parameters 
in the optical spectral range. In this paper we are for the 
first time in position to determine galaxy kinematics us- 
ing SSP templates in the near-IR Call spectral range. It 
is generally assumed, that kinematic stellar templates are 
required to be obtained with the same ins trumental set- 
ting as the one used for the galaxy spectra ( Sargent et al.l 
Il977t iTonrv fc Da,vislll979h . In this paper we show that 
this is no longer a requirement, at least at the resolutions 
of the spectra considered here. 

This paper intends to be a first step in the exploration 
of SSP models as templates by comparing velocity pro- 
files, for a sample of 19 galaxies, with those obtained with 
our best stellar or mixture of stellar templates. We show 
that in all cases SSP models produce kinematics of at 
least at the same level of quality, and are quite often 
better than the best stellar templates. In order to trace 
correlations between kinematic and population features 
along the galaxy mino r axes, we measure the CaT* index 
(ICena,rro et, a1.ll2001ah. 

Sect. 2 summarizes the observations and data reduc- 
tion. In Sect. 3 we describe the different types of template 
samples and the way the kinematical profiles have been de- 
termined, including a detailed comparison between stellar 
templates and SSP models. The obtained kinematics is 
discussed in Sect. 4. Kinematical profiles are graphically 
presented in Appendix EI Tables containing the numer- 
ical values presented in appendix E] will be available in 
electronic form. 



2. Observations & Data Reduction 

2.1. Sample 

The sample comp rises 20 galaxies of types SO-Sbc se- 
lected out of the iBalcells Peletie'il (|l994l) sample of 
bulges. The latter is a diameter-limited sample of galaxies 
with inclinat ions above 50 degree from the UGC catalog 
(|Nilsonl[l9 95): it is has been the subject of se veral stud- 
ies on bulges (see IFalcon- Bar roso et alJ 120021 and refer- 
ences therein). The 20 objects were those we had observed 
with HST/WFPC2 and NICMOS, hence high-spatial res- 
olution images are available for them at B^I^H bands 
([Peletier et al.lll99fJ) . The main properties of the sample 
are presented in Table [l] 



2.2. Observations 

We obtained long-slit spectra for our sample of galaxies at 
the 4.2m William Herschel Telescope of the Observatorio 
del Roque de los Muchachos at La Palma between 11 and 
13 July 1997. Twenty galaxies were observed between 8360 
- 9170 A using the ISIS double spectrograph. We discarded 
NGC 5577 as the spectra did not have sufficient signal-to- 
noise. We are thus left with 19 objects (see Table GJ. The 
red arm was equipped with the Tektronix (1024x1024) 
TEK2 CCD (0.36"per 24 fim pixel) and the R600R grat- 
ing, providing a spectral resolution (FWHM of arc lines) 
of 1.74 A (59.6 km/s) with a spectral dispersion of 0.8 
A/pix. We used the 6100 A dichroic for simultaneous red 
and blue arm exposures. In this paper we discuss the red 
arm spectra only. Typical exposure times per galaxy were 
1500 sec. The length of the slit was 3.6'. The slit width was 
1.2", matching the seeing at the time of the observations. 
Arc line exposures were taken before and after each tar- 
get exposure. Tungsten continuum lamp exposures were 
taken with the red arm after each target exposure, for 
fringe calibration. Twilight sky exposures were taken ev- 
ery night for flat fielding. We also obtained spectra for 11 
stella r templates of types A 2 to M4, from the Lick list of 
stars (|Worthev et al.lll994h . listed in Table 

2.3. Data Reduction 

The data reduction was carried out using the IRAF 1 pack- 
age. The general standard procedures were applied to the 
data (bias subtraction, flat fielding, fringe correction, sky 
su btraction, wavelength and flu x calibration), as described 
in IFalcon- Bar roso et al.l ((2002). Cosmi c ray removal was 
performed using the REDmE package ()Cardiell ll999). 

3. Data analysis 

3.1. Stellar Templates 

In Table [21 we list the observed stellar templates. A repre- 
sentative set of spectra is shown in Fig.Q] The wide range 
of spectral types covered allows us for an exploration of 
the best templates for matching galaxy spectra. The most 
apparent variation from top to bottom in Fig. [l] is the de- 
crease of the Paschen lines toward later types. P14, in a 
region devoid of other features, is a strong discriminant of 
the presence of stars earlier than G. P13 and P15, when 
present, give characteristic asymm etric profiles to the sec - 
ond and third lines of the Call ([Cenarro et al.l l2001bh . 
Following previous wo rks in the literature (jRix fcW hite 
1992; lHau et all 1999) we have included in our analy- 
sis a mix-template from a linear combination of stars in 
our sample. However we have used a different approach 
to derive our optimal mix- template. Given the impor- 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 
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Table 1. Galaxy Sample 



NGC 


V LG 


Scale 


Type 






B/D 


PA 


Exp. Time 


6 D 


e B 


^bulge 


Best Tempi. 


S/N 
pix -1 




km/s 


kpc/" 




mag 


mag 




deg 


sec 






II 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ii) 


(12) 


(13) 


(14) 


5326 


2576 


0.25 


1 


-22.16 


-21.22 


0.73 


-138 


1500 


0.50 


0.43 


5.97 


+0.2, 3.16 


83 


5389 


1990 


0.19 





-21.32 


-20.70 


1.30 


93 


1500 


0.80 


0.46 


2.85 


Ml III / +0.2, 4.47 


61 


5422 


1929 


0.19 


-2 


-21.96 


-21.33 


1.28 


64 


1500 


0.80 


0.35 


3.33 


+0.2, 7.08 


68 


5443 


2060 


0.20 


3 


-21.53 


-19.09 


0.12 


-52 


1400 


0.68 


0.54 


1.98 


+0.2, 3.16 


37 


5475 


1861 


0.18 





-20.98 


-18.72 


0.14 


-101 


1500 


0.68 


0.24 


2.45 


+0.2, 1.41 


57 


5587 


2291 


0.22 





-21.07 


-19.43 


0.29 


73 


1500 


0.70 


0.47 


1.51 


-0.4, 14.12 


36 


5689 


2290 


0.22 





-22.29 


-21.47 


0.89 


-5 


1500 


0.75 


0.48 


2.04 


0.0, 4.47 


58 


5707 


2354 


0.23 


2 


-21.43 


-20.11 


0.42 


-59 


1500 


0.75 


0.20 


3.79 


0.0, 2.51 


51 


5719 


1684 


0.16 


2 


-21.60 


-20.00 


0.30 


5 


1500 


0.64 


0.61 


7.71 


K3 lb / -0.4, 3.55 


51 


5746 


1677 


0.16 


3 


-22.68 


-21.70 


0.68 


-99 


1200 


0.84 


0.52 


2.37 


Ml III / +0.2, 3.16 


44 


5838 


1337 


0.13 


-3 


-21.89 


-20.93 


0.71 


134 


1500 


0.65 


0.22 


8.39 


0.0, 12.59 


86 


5854 


1708 


0.17 


-1 


-21.43 


-20.18 


0.46 


-33 


1500 


0.70 


0.32 


4.85 


-0.7, 6.31 


64 


5879 


1065 


0.10 


4 


-20.30 


-19.08 


0.48 


89 


1500 


0.62 


0.32 


0.85 


-0.4, 3.16 


41 


5965 


3603 


0.35 


3 


-22.92 


-21.76 


0.53 


-38 


1500 


0.84 


0.50 


2.32 


-0.4, 5.01 


46 


6010 


1923 


0.19 





-21.57 


-19.88 


0.27 


12 


1500 


0.75 


0.28 


2.70 


0.0, 14.12 


67 


6504 


4680 


0.46 


2 


-24.57 


-23.30 


0.45 


135 


1500 


0.83 


0.44 


4.05 


-0.4, 3.16 


52 


7331 


1138 


0.11 


3 






0.001 1 


80 


1500 


0.63 1 


0.62 


4.57 


-0.4, 15.85 


117 


7332 


1550 


0.15 


-2 


-21.86 


-20.52 


0.41 


-114 


1500 


0.74 


0.29 


2.99 


-0.4, 11.22 


101 


7457 


1114 


0.11 


-3 


-20.91 


-20.91 


10.00 


35 


1500 


0.48 


0.34 


3.29 


+0.2, 2.51 


51 



Notes: Description of the columns: 

(1) : the NGC numbers of the galaxies. 

(2) : Recession velocity of each galaxy in km/s, corrected to the Local Group (iKarachentsev & Ma karovl 1 1996^ 1 . from optical 
heliocentric velocities given i n the Third Reference Cata logue of Bright Galaxies (|l99 J ) (hereafter RC3). 

(3) : Scale in kpc/arcsec from iFalcon- Bar roso et all (120021) . 

(4) : Type index T from the RC3. 

(5) and (6): Absolute magnitudes, in i?-band, from IPeletier Balcellsl (|l997l) (Hp= 50 km /s/Mpc). 

(7) : Bulge-to-disk ratio from an i?-band bulge-disk decomposition follo wing Kent's (I1984T) metho d (jPeletier Balceilslll997l) . 

(8) : The position angle, in degrees, (N - E) of the dust-free minor axis (Pele tier Sz Balcellslll997F) . 

(9) : Integration time, in s ec, for each galaxy during the observations. 

(10) : Disk ellipticity from IPeletier fc Balcellsl (I1997T) . 

^(11 ^ Effective bulge ellipticity, from HST i7-band ellipticity profiles, measured at the geometric effective radius ( Balcells et alJ 

(12) : Minor axis bulge radius, derived from iBalcells et alJ (|2003f) HST if-band profiles, as the radius at which the contri- 
bution from the disk an d bulge light are equal, except for NGC 7332 for which a ground-based image was used instead 
(Pele tier Balceilslll997l) (Alister W. Graham, private communication). 

(13) Best template fit for each galaxy. The metallicity ([m/H]) and age (Gyr) of the best fit SSP models is shown. In those cases 
where the stellar template performed better, the adopted SSP best template is also given. 

(14) Sign al-to-noise ratio p er pixel in the central aperture (1.2 x 0.36 arcsec 2 ). 
1 From lPrada et all (Il996l) 



tance of the continuu m shape in this wavelength region 
(ICenarro et al .l l2001at IVazdekis et ail fcoO.U we have cre- 
ated an optimal template that provides the best fit to the 
overall continuum shape of each galaxy (i.e. fitting the re- 
gions outside the main Call features, but including the 
TiO bands), before broadening with the appropriate kine- 
matics. It is important to perform this fit using a com- 
bination of unconvolved stellar spectra, since any input 
kinematics comes from a particular template and could 
bias our result. In most cases the optimal mix-template is 
formed by G8III, K3Ib and M4V stars. The importance of 
the M4V star in the final template is justified by its abil- 
ity to match the TiO molecular bands. As we will show 
in Sect. 13.41 this mixture of stellar templates helps to al- 



leviate the template mismatch problem and it is an ideal 
benchmark to test our SSP models. 

3.2. SSP model templates 

We have used 78 SSP synthetic models at FWHM=1.5 A 
(0.85 A/pix) in the Call triplet region, corresponding 
to a single-age, single-metallicity stellar populations of 
[Fe/H]=-0. 7,-0.4, 0.0, and +0.2. For each metallicity, 
the age varies from 1 to 17 Gyr, except for [Fe/H]=— 0.7 
starting from 6 Gyr and for [Fe/H]=— 0.4 starting 3 Gyr. 
These limitations arise from the limited coverage of atmo- 
spheric parameters of the stellar sample used by the m od- 
els. More details can be found in lVazdekis et al.l (|2003l) . In 
order to adapt these models to the instrumental resolution 
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Table 2. Stellar Templates 



Name 
(1) 


a 


(h m s) 
(2) 


S (d m s) Spectral Type M v 
(3) (4) (5) 


BARNARD 1 


17 


57 


48.50 


+04 


41 


36.2 


M4 Ve 


9.54 


GL818 2 


21 


05 


19.75 


+07 


04 


09.5 


K6 V 


8.33 


HR0072 


00 


18 


41.87 


—08 


03 


10.8 


GO V 


6.46 


HR4521 


11 


46 


55.62 


+55 


37 


41.5 


K3 III 


5.26 


HR5340 


14 


15 


39.67 


+19 


10 


56.7 


K3 lb 


1.23 


HR5681 


15 


15 


30.16 


+33 


18 


53.4 


G8 III 


3.47 


HR5826 


15 


31 


24.93 


+77 


20 


57.7 


K5 III 


5.00 


HR5854 


15 


44 


16.07 


+06 


25 


32.3 


K2 III 


2.64 


HR6685 


17 


55 


25.19 


+26 


03 


00.0 


F2 Ibe 


5.47 


HR8334 


21 


45 


26.93 


+61 


07 


14.9 


A2 lab 


4.31 


HR8795 


23 


07 


00.26 


+09 


24 


34.2 


Ml III 


4.55 



Notes: 

(1) HR names from the Bright Star Catalogue 
ftToffleit WarrenlllQQ^. 

(2) , (3) J2000 coordina tes from SIMBA D 

(4) , spectral type from Worthev et a Q (1 19941) 

(5) from SIMBAD (http://simbad.u-strasbg.fr/Simbad) 

1 NS V 9910 from the New catal ogue of suspected variable 
stars (iKukarkin fc Kholopovlfl982h 

2 HP 200779 from the Henr y Draper Catalogue and Extension 
(|Cannon & Pickerin j|jaa3) 



P18P17P16 P15 P14 P13 P12 A2Iab (HR8334) 




8400 8500 8600 8700 8800 8900 

A [A] 

Fig. 1. Representative sample of continuum subtracted 
stellar templates used in this paper. The most important 
features found in those spectra have been indicated: the 
Paschen series (from Pll to P18), the Ca n triplet (Cal, 
Ca2 and Ca3) and several metal lines (Fel, Mgl and Til). 



of our data, we have convolved them with a Gaussian of 
<j=0.38 A. This procedure relies on the assumption that 
a Gaussian represents a good match to the instrumen- 
tal PSF (Point Spread Function) of the observations. We 
have tested whether the arclines followed a gaussian shape, 
finding that indeed they do, with only small deviations 
(b elow 0.01) in the h% an d fe 4 p arameters, as described 
bv lvan der Marel fc Franxl f|l993h . If J13 and h± had been 
found to be significantly different from 0, we would have 



had to convolve the models with the line profile of the 
arc-lines. 

3.3. Extraction of the kinematics 

Kinematic profiles were derived from the data in the re- 
gion of the near-IR Call triplet (8498 A, 8542 A, 8662 A). 
For each galaxy, we generated two sets of coadded galaxy 
spectra, with minimum signal-to-noise ratios (S/N) of 33 
and 20 per pixel, respectiv ely. We made use of the program 
FOURFIT developed bv Ivan der Marel fc FramJ (|l993), 
available from the first author's website. A Gauss-Hermite 
fit to the LOSVD was performed to obtain mean radial 
velocities, velocity dispersions, hs and /14 moments for 
the data with signal-to- noise ratio above 33 per pix -1 . 
For data with signal-to-noise ratio between 33 and 20 per 
pix -1 a pure gaussian fit was used. The choice of S/N is 
explained in Sect. 13.51 

3.4. Selection of the best fitting template 

We have determined the kinematic profiles for our sample 
using each of the 90 templates we had available (78 SSP 
models, 11 stellar templates and the best mix-template 
for each galaxy). We have considered as our best fit 
template the one that minim ized the x 2 7 as defined in 
Ivan der Marel fc Franxl ([1293), for the central position of 
the galaxy profile, which has a S/N well above 33. One 
of the features of this code is the possibility of filtering 
the input spectra in Fourier space so that any remain- 
ing low-order component not removed by the continuum 
subtraction or high-frequency noise can be filtered out. 
The resulting x 2 is then evaluated over a finite range of 
wavenumbers determined by the filtering parameters. For 
our dataset the results from FOURFIT were weakly de- 
pendent on the lower cutoff and insensitive to the choice 
of the upper cutoff (given that we chose the same upper 
cutoff for all the galaxies and that this cut was well above 
the sigma of the galaxy). We also reduce the uncertain- 
ties in the derived kinematics by flux-calibrating the data 
and the template stars before applying the continuum sub- 
traction because the normalizing polynomia for the data 
should be very similar to the ones used for the templates, 
otherwise, the normalization procedure might introduce 
non-negligible variations in the spectra, including in the 
equivalent widths of the lines. 

We establish whether SSP, stellar templates or mix- 
template provide better fits to the galaxy spectra by com- 
paring the global reduced x 2 of the fits. A typical distribu- 
tion of x 2 values is shown in Fig. [21 The lines trace reduced 
X 2 for the SSP models, labeled by age and metallicity, 
while the crosses are the reduced x 2 of the stellar tem- 
plates labeled by their spectral types and combined stellar 
template. The patterns in this Fig. are representative of 
those found for all the galaxies in the sample. There is a 
group of stellar templates that do significantly better than 
the others; in particular, the A2, F2 and M4 templates 
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5 1015 A2I F2I GOV G8III K2II K3I K3II K5II K6V M1III M4V Mix 
Age [Gyr] / Stellar Templates 

Fig. 2. x 2 distribution for NGC 5326. Open square indi- 
cates the minimum x 2 for the best SSP fitting model, from 
which kinematics has been derived for this galaxy, x 2 val- 
ues for the stellar templates A2 lab, F2 Ibe, GO V and M4 
Ve have been plotted with arrows indicating that they fall 
outside the plot limits. Their x 2 values are: 2.197, 1.313, 
0.613, 1.715 respectively. 



consistently yield poor fits in all the galaxies. For the SSP 
models, the x 2 distributions are fairly flat as a function of 
age for a given metallicity, except for ages below ~2 Gyr, 
where the quality of the fits quickly degrades. This behav- 
ior is expected since, at younger ages, the Call equivalent 
widths significantly decrease and the Pas chen series be- 
come noticeable (see IVazdekis et all 120 03). At each age, 
the fits are sensitive to metallicity. 

For the spectrum used in Fig. [2J the best fit is ob- 
tained with a SSP model corresponding to [Fe/H]=+0.2 
and 3.16 Gyr. It is worth noting that this SSP model 
should not necessarily be considered as the best fit to the 
stellar populations of this galaxy (see below). For the en- 
tire sample, the x 2 f° r the best SSP model, best stellar 
template and mix-template are quite similar, and for 16 
of the 19 objects the SSP templates provide marginally 
better fits. This result demonstrates that synthetic SSP 
spectra, convolved to the resolution of the data, can pro- 
vide accurate templates for kinematic analysis of galaxy 
spectra. The high quality of the SSP fits is most-likely 
due to their ability to minimize the template mismatch 
problem. While the stellar populations of bulges are likely 
to contain a mixture of ages and a range of metallici- 
ties, the SSP models provide a better approximation to 
galaxy spectra than a single stellar spectrum. Our results 
indicate that late Kill stars and early M stars provide 
a better approximation to the integrated light of bulges 
than A , F, and G stars in agreement with IVazdekis et al.l 
(|2003l) . Tabled Col. 13, lists the overall best template 
for each galaxy. Note that, for three galaxies (NGC 5389, 
NGC 5719, NGC 5746), the M1III and K3I stellar tem- 
plates provide the overall best fit to the galaxy spectra. 

In Fig. 03 we plot the central spectrum of NGC 5326 and 
the best-fitting stellar (Ml III), mix and SSP templates, 
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Fig. 4. Comparison of the profiles obtained with our 
best stellar templates (G8III, K2III, K3I, K3III, K5III, 
K6V, M1III), mix-template and best SSP model fit of 
NGC 5326. 



convolved with their respective kinematic solutions. The 
SSP model represents a very good fit for the whole sam- 
ple of galaxies spectrum. Fitted output parameters for the 
best template stars (G8III, K2III, K3I, K3III, K5III, K6V, 
M1III), mix-template and best SSP model are shown in 
fig.^Ifor the entire profile of NGC 5326. The results show 
an excellent agreement. Therefore we adopt the solution 
from our best SSP template for each galaxy to obtain kine- 
matic profiles. 

A point to remark is that no inferences on the stellar 
population age and metallicity of the galaxies are pre- 
sented here from the best-fitting SSP models. The reason 
is that an important relation has been shown to exist be- 
tween the continuum shape and the stella r populations in 
this spectral region dCenarro et, a ll l2003t IVaz dekis et, all 
2003). Since we remove this continuum, as required by 
FOURFIT, we neglect its valuable information for con- 
straining galaxy stellar populations. Therefore the most 
useful information left in our spectra is the equivalent 
width of the Call feature, which does not vary much for 
SSPs of the age and metallicity regimes considered here. 
The wavelength region that we have used to derive the 
kinematics excludes the strong TiO molecular bands seen 
in this spectral range. Therefore we are left with an almost 
featureless continuum shape which is straightforward to 
subtract. A different methodology is required to be able to 
analyze their stellar populations on the basis of the SSP 
models. Methods for extracting use ful information from 
the full SSP spectra are di scussed in IVazdekis! (|l999h and 
IVazdekis fc Arimotol f|l999h . The stellar population analy- 
sis of t his galaxy sample will be pre sented in a forthcoming 
paper ((Falcon- Bar roso et aLlEomh. 

The use of SSP models has a potential to minimize the 
template mismatch problems that have traditionally af- 
fected galaxy kinematic measurements. Whereas the Kill 
stellar templates have been shown to provide reasonably 
good fits to elliptical galaxy spectra in the optical region, 
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Fig. 3. Comparison between our best stellar template, mix-template, and best SSP model fits, for the continuum 
subtracted and low-order filtered, binned central spectrum of NGC 5326. The fitting region has been indicated with 
a line at the bottom of the Fig. In this case the best fitting template is the SSP model. 



one expects that the SSP models should significantly im- 
prove the kinematical fits. The blue and visible have been 
proven to have large sensitivi ties to the age, metallicity 
and abundance ratios iRose et al.l Il994t IWorthev et al.l 
Il994t IVazdekis et "all 1 19961: IVazdekieJ 11999^ . Our results 



suggest that the requirement of observing template stars 
with the same instrumental setting as the galaxies is not 
compelling. The use of SSP models as templates allows 
us to use all of the available observing time for the tar- 
gets, without the need to take spectra of template stars. 
An important application of SSP models should be in the 
field of high-redshift galaxy kinematics, where the mea- 
surement of the same spectral range in stars and target 
galaxies would require changing the instrumental setup, 
and where the redshift of the target objects is not gener- 
ally known a priori. 



The main limitation of the use of SSP models is that 
it relies on a well designed stellar library to produce the 
full stellar population synthesis. However, the method is 
expandable to include combinations of SSP models to ac- 
count for composite populations. Potential applications 



are the study of age differences between bulge and disk, or 
population signatures of kinematically-decoupled cores. 



3.5. Errors & Simulations 

Although FOURFIT gives the formal errors of the fit, we 
have performed simulations to establish the level of con- 
fidence the code can produce on (V, a, hs and h±) at 
S/N=33 pix- 1 , and (V, a) at S/N=20 pix" 1 . 

For the simulations two different tests have been per- 
formed. In the first one we convolved an SSP model 
([Fe/H]= 0.0, t= 5.01 Gyr) with a wide range of LOSVDs 
to generate model galaxy spectra. The LOSVDs included 
several velocities, velocity dispersions, I13 and I14 mapping 
a representative parameter space covered by our sample 
of galaxies (V E [800,2000], a e [50,300], h 3 and h 4 G 
[-0.1,0.1]). The inclusion of higher velocities (up to 4000 
km/s, to match our furthest galaxy) doesn't affect our re- 
sults. Eight different samples of white noise were added to 
our model galaxy spectra to yield S/N ratios of 20 and 33 
pix -1 . For each one of the noise sample at the two S/N 
levels the full battery of LOSVDs was recovered. Given 
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Table 3. Errors from simulations 



Par am. 


V [km s- 1 ] 


a [km s l ] 


h 3 


h 4 


S/N 


20 33 


20 33 


20 33 


20 33 


std(SSP) 


15 13 


29 15 


- 0.07 


- 0.08 


std( Stars) 


18 16 


34 18 


- 0.07 


- 0.10 



Notes: Standard deviations for each parameter at the differ- 
ent signal-to-noise ratios measured from the simulations (see 

Sect. 023- 



the fact that all SSP models give very similar x 2 values, 
only a representative subsample of 3 SSP models covering 
the models parameter space ([Fe/H]=-0.7, Age=7.98 Gyr; 
[Fe/H]=0.0, Age= 5 Gyr; [Fe/H]=0.2, Age=15.85 Gyr) 
was used to perform these simulations. We have simu- 
lated both Gaussian and non- Gaussian LOSVDs, however 
we only fitted pure gaussians to the spectra of S/N=20 
pix -1 , for which we only modeled the V and a errors due 
to LOSVD mismatching in the S/N=20 spectra. In the 
second test we created the galaxy spectra from a stellar 
template (K3I). Spectra from the K3I, K5III and M1III 
stars were used as kinematic templates. The choice of 
these templates is due to their ability to produce, overall, 
the lowest values of x 2 when fitted to the sample of galax- 
ies. We have convolved the model galaxy spectra with the 
same LOSVDs and added the same white noise samples 
as in the first test. 

The main results of the simulations are summarised 
in Table 01 The Table gives the standard deviations of 
the distributions (parameter^ - parameter^), obtained 
from the simulations in each test, for the full set of pa- 
rameters and white noise samples. This is a measurement 
of the uncertainty when working at the mentioned S/N 
ratios. We use the uncertainties obtained from the sim- 
ulations with the SSP models and display them as error 
bars on the right-hand side of each kinematic profile, see 
Appendix A. In Fig. the left panels show the distri- 
bution of the recovered LOSVDs compared to the input 
ones at each sampled value of a and I14 for each test for 
S/N>33. The values shown below each histogram repre- 
sent the median (within square brackets) and standard 
deviation of the distribution, found at those values of a 
and h4. The right panels, on the other hand, show the 
total distribution and standard deviation (i.e. the sum of 
the individual ones) from which we have obtain the values 
shown in Table 01 No significant differences are found be- 
tween the distributions obtained from the SSP models and 
stellar templates. This result strengthens our approach of 
using SSP models as kinematic templates. 

The simulations show two effects worth mentioning. 
First, there is no systematic underestimate of the velocity 
dispersion for low input velocity dispersions (see Fig. EJ). 
This effect, if present, might bias our velocity dispersion 
results for our two smallest bulges, s ee Fig. El Such an ef- 
fect has been previously reported bv I Joseph et al.l ()200ll) 
using a different code (i.e. FCQ), finding an overestimation 
of ^10 km/s at a=50 km/s. Second, as we move towards 
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Fig. 5. Recovery of a and h 4 via FOURFIT using SSP 
models and stellar templates. The left panels show the 
distribution of the recovered LOSVDs compared to the 
input ones at each sampled value of a and I14 for each test 
at S/N>33. The right panels, on the other hand, show the 
total distribution (i.e. the sum of the individual ones on 
the right panels). For each a and I14 we plot a histogram 
of the A(cr)=cr(out) — cr(in) and A(h4)=li4(out) — li4(in) 
obtained from the 14400 realisations corresponding to 8 
white noise realizations. The values shown below each his- 
togram represent the median (within square brackets) and 
standard deviation found at those values of a and I14. 

larger velocity dispersions, the uncertainty in <j also in- 
creases (from 6 km/s at <r=50 km/s, to 15 km/s between 
a=100-200 km/s and 32 km/s between a=250-300 km/s). 
Hence, although the Gauss-Hermite expansion conforms a 
description of the LOSVD with a set of orthogonal func- 
tions and therefore we only expect major coupling betw een 
7 (the peak of the LOSVD) and a (jLarsen et alJll983h . V 
and /13, and a and /14, we have detected other correla- 
tions when working at low signal-to-noise ratios. The ex- 
pected coupling of V with ft 3 and a with appears in our 
simulations. Similar results have been found by authors 
who have performed s imilar tests ((van der Marel fc~F ranx 
ll993tlHau et al.lll999h . 

We have decided to impose a S/N=33 as the min- 
imum level to determine ft3 an d /14. This is a slightly 
lower level than that suggested bv lvan der Marel fc Franxl 
(1993) (S/N-40 pix" 1 ) to give a good estimation of 
the Gauss-Hermite higher order moments. Nevertheless 
in most cases we obtain S/N>40 for the profile central 
position (see Table GJ. 

3.6. Measurement of the Ca~H index 

We ha ve measured the rece ntly defined Paschen-free CaT* 
index ((Cenarro et al.ll2001at) . that measures the strength 
of the Call lines corrected from the contamination by 
Paschen lines, from our flux calibrated spectra. We have 
made used of the FORTRAN routine indexf its (part of 
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Fig. 6. Comparison betw een our best SSP model fit and 
IVega Beltran et, al.1 fcOOlh for NGC 5854. 



the REDmE package), kindly made available to us by 
Javier Cenarro. No velocity dispersion correction and zero- 
point offset have been performed. Errors in the CaT* in- 
dex include contributions from radial velocity uncertain- 
ties and from S/N. The first one is obtained by the CaT* 
index from 100 random simulations within the velocity 
range [V — AV,V -\- AV]. The latt er is determined using 
eq. A37 from lCenarro et ah! ((2001 al) . Both terms are added 
in quadrature. The uncertainty from the S/N dominates 
the CaT* error. 



Call triplet index CaT*, along the galaxy's minor axis. 
Abscissae are centered on the brightness peak of the long- 
slit spectrum, i.e. on the photometric center of the galaxy 
since we don't expect large centering errors. Negative ab- 
scissae correspond to the minor axis that projects away 
from the disk (the "dust-free" side), whose position an- 
gle on the sky is given at the top next to the galaxy's 
name. Circles are measurements derived from S/N > 33 
pix -1 spectra, and triangles are V, a and CaT* derived 
from 20 < S/N < 33 spectra. Error bars in each point are 
the formal error estimates coming from FOURFIT code. 
Measurement uncertainties derived from the simulations 
described in Sect. 13.51 are shown with vertical bars on the 
right edge of each plot for the S/N > 33 spectra and, for 
V and a, for 20 < S/N < 33 as well (longest bar). The 
contour plot draws isophote s of the HST NICMOS F160W 
images ((Peletier et all l999). drawn to scale. The length of 
the 1.2 arcsec wide slit delimits the range along the minor 
axis profile from which spectra has been binned to obtain 
a S/N>20. The high spatial resolution of the NICMOS im- 
ages (0.075 arcsec/pixel, PSF=0.18 arcsec FWHM) pro- 
vides a clue on the level of substructure existing within 
the spectroscopic slit. Due to the lack of a NICMOS image 
for NGC 7332, a WFPC1 image, retrieved from the HST 
archive, has been used instead. Vertical dotted lines indi- 
cate the minor axis bulge radius and center of the galaxy 
(see Table [I] for more details). 

Owing to the short integration times and the S/N re- 
quirements for kinematic extractions derived from the sim- 
ulations (Sect. 13.5)) . the kinematic profiles typically ex- 
tend to 5-10 arcsec from the nucleus only. Only for the 
large galaxy NGC 7331 do kinematic profiles extend over 
20 arcsec on the dust-free side. All profiles are drawn to 
the same horizontal scale. In the following subsections we 
briefly address the global behavior of each kinematic pa- 
rameter. 



3.7. Comparison with the literature 

We have compared our best SSP model fits with those in 
the literature. Unfortunately, published minor axis kine- 
matic profiles are scarce. We present the comparison in 
Fig. and □ for NGC 5854 and NGC 7332 respectively. 
An excellent agreement is found for the two galaxies, with 
no major systematic deviations, supporting our approach 
of using synthetic SSP models. This agreement is par- 
ticularly good for the velocity dispersion of NGC 5854, 
where our results follow the same trend obtained by 
IVega Beltran et all ()200lh . The c omparison for NGC 7332 
reveals also good agreeme nt with iFisher et al.l (|l994l) and 
ISimien fc PruiHc] Eli3) • 

4. Results 

Kinematic profiles for the 19 galaxies in the sample are 
shown in Fig. A.l (Appendix A). For each galaxy we 
show profiles of mean heliocentric velocity V, velocity dis- 
persion (a), hs and J14 Gauss-Hermite coefficients, and 



4.1. Velocity profiles 

As expected for the minor- axis alignment of the slits, mean 
velocity profiles are fairly flat. However, the profiles are 
not entirely featureless. We note three types of features: 

4.1.1. Rotation in the outer parts 

We detect non-zero rotation in the outer parts of the pro- 
files (opposite mean velocities at both sides of the center, 
at the 2-cr level) in 6 out of the 19 galaxies (NGC 5389, 
5689, 5746, 5838, 5854, 7457). Typical rotation amplitudes 
are in the range 10 < V < 20 kins" 1 . Such estimates of 
velocity amplitude and frequency of occurrence are highly 
dependent on the depth of the data. For two galaxies 
(NGC 5475, 5587) the profiles do not reach 5 arcsec and 
contain no information on rotation at large radii. Brighter 
galaxies in which the spectra extend farthest from the cen- 
ter show the highest V at large radii. In most of these 
cases, the outer velocity measurements sample disk or 
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Fig. 7. Comparison between o ur best SSP model f i t with 
iFisher et all ((l994h (F94) and LSimien fc Prugniel (|l997l) 
(SP97) for NGC 7332. 

bulge-disk transition regions. The non-zero velocities indi- 
cate that our spectroscopic slits were not perfectly aligned 
with the line of nodes. The slit was placed orthogonal to 
the major axis of the outer isophotes. Hence the misalign- 
ment might reflect an isophote twist in the outer parts 
sampled by the spectra, due eg. to warps or spiral arms, 
or a distortion of the inner disk velocity field due to a bar 
or a spiral arm. Note that a misalignment of 5° or more 
could easily explain such amplitudes in velocity. 

4.1.2. Inner rotation 

We detect inner minor-axis rotation in 8 galaxies 
(NGC 5443, 5587, 5689, 5707, 5838, 5965, 6010, 7457) and 
possibly also in NGC 5389, 5422, 5854. Whenever the pro- 
files extend beyond the nuclear region, inner rotation is 
delimited by a change in the velocity profile slope. These 
bulges thus contain kinematically distinct regions (here- 
after KDR). We do not intend to set a link between these 
KDR and counterrotating or otherwise kinematically de- 
coupled cores in ellipticals and bulges; full counterrotation 
cannot be established from our minor-axis kinematic pro- 
files. Rather, we associate the velocity features to slight 
variations in the kinematic structure, related to inner com- 
ponents such as disks or bars. 

Theoretically the simplest way to explain a minor- 
axis rotation in a bulge is to assume we are dealing 
with a triaxial object. In a triaxial object the position 
angle of the apparent minor axis can be different from 
that of the projected in trinsic short axis of the object 
(|de Zeeuw fc Franxlll99lh . Moreover, the photometric and 



kinematic minor axes may be misaligned as a result of 
the wide range of possible locations of the total angu- 
lar momentum in a triaxial object (Bin nevl Il985l) . In 
practice, whether minor-axis rotation traces a triaxial 
bulge can only be che cked using integra l field spectro- 
graphs e.g. SAURON (iBacon et all BlOn. or measuring 
kinematics along different posi tions angles of the galaxy 
(jfttatler fr Pmecker-HanejEjaj). 

The typical minor-axis radial extent of the KDR is 
2-3 arcsec, or 0.5 kpc at the median distance for the 
sample. In all of the eight galaxies listed above as KDR, 
isophotal analysis of the HST/NICMOS F160W images 
(our work, unpublished) reveals an inflection in the C4 
isophotal coefficient of the cos(4#) term that measures 
disky deviations of th e isophotes from pure ellipses (e.g. 

recense n et al.l ll992 N ). C4 amplitudes of the features are 
above 1% in all cases. This suggests the presence of flat- 
tened subcomponents in the inner regions of the bulges 
showing up as excesses in t he bulge surface brightness 
profiles ((Balcells et al.l 120031) . One probably sees a mix- 
ture of inner disks and inner bars. In four of the eight 
cases the isophotal analysis reveals changes in isophote 
position angle (APA > 10°) at the same radius, which 
are likely to correspond to inner bars but could also trace 
misaligned disks. The small rotation amplitudes and the 
lack of isophotes elongated perpendicularly to the disks 
suggests that none of our bulges show the t ype of orthog- 
onal bulge kinematics found in NGC 4 698 (iBertola et all 
Il999l) and NGC 4672 (jSarzi et al.ll2000h . The four remain- 
ing cases with no isophotal misalignment probably corre- 
spond to inner bars, as an aligned disk has its line of nodes 
parallel to the minor axis. 

The most prominent nuclear KDR is that of 
NGC 5838. The nucleus of NGC 583 8 harbors an ordere d 
ring or disk of dust and young stars ((Peletier et al.l ll999). 
which is oriented orthogonal to the slit. The dust may be 
responsible for the asymmetric rotation pattern. 

It is worth noting th at none of the KDR have distinct 
colors B — I and I — H ((Peletier et al.lll99^ on the dust- 
free side of the galaxies, nor do they have corresponding 
features in the CaT* profiles (Fig. Al). This may suggests 
that the velocity featu res are not associa ted to distinct 
stellar populations (e.g. lDavies et al1l200ll) . 

4.1.3. Kinematically offset nuclei 

In two galaxies (NGC 5707, 7331) the velocities at large 
radii appear offset in the same direction relative to the 
nuclear velocity. NGC 5707 has a nuclear disk or bar ex- 
tending over 1 arcsec, which is well aligned with the ma- 
jor axis of the outer disk and partially obscured by a dust 
lane; the dust might be responsible for the offset velocity 
measurement at the photometric center of the galaxy. In 
NGC 7331, the entire region from -4 arcsec to +4 arcsec 
appears displaced 15-20 krns" 1 from the velocities further 
out. It is unclear whether dust extinction can be responsi- 
ble for this structure, or whether the bulge might actually 
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wander in the potential of the galaxy. The major-axis kine- 
matics of this galaxy is complex, and its outer spiral struc- 
ture is asymmetric; this could indicate a recent accretion 
event and/or the motion of the lumin ous material within 
the d ark matter halo of the galaxy ([Noordermeer et "all 
Results on th e major- axis kinematics of this galaxy 
are puzzling. While IPrada et all (^996) reported to have 
found the bulg e in NGC 7331 to counter-rotate, a more 
recent work bv iBottemal (|l999l) contradicts this result. 



4.2. Velocity dispersion profiles 

Velocity dispersion profiles trace the dynamical tempera- 
ture of the bulges and are key ingredients for dynamical 
modeling of the bulges. They also provide diagnostics for 
distinguishing between bulge and disk on the basis of the 
internal dynamics. 

For our galaxies, the minor-axis velocity dispersion 
profiles show a wide diversity of shapes, from flat to 
strongly centrally peaked. Velocity dispersion gradients, 
(Act per unit minor-axis r e ) were computed by fitting lin- 
ear relations to the negative-abscissa dispersion profiles 
of Fig. A.l, which correspond to the dust-free side of the 
bulges. Minor- axis r f were derived using the r e values from 
iBalcells et alJ ([20031) . which derive from i^-band surface 
brightness profiles from i^5T-NICMOS images, and ellip- 
ticities at the bulge effective radius €B{r e ) obtained from 
the ellipticity profiles in the same paper. Dispersion gra- 
dients range from to -70 km/s/r e , with a median of -21 
km/s/r e and a mode of only ~10 km/s/r e . 

Fig. [HI displays all the velocity dispersion profiles, pre- 
sented on two panels for clarity; the top/bottom panel 
shows the half sample with highest /lowest V(<r). In or- 
der to show whether the profiles sample only the bulge or 
extend to the disk region, we normalize each abscissa to 
f 'bulge, the minor-axis radius where bulge and disk con- 
tribute equally to the surface brightness profile (listed 
in Col. 12 of Table Most profiles extend well beyond 
r/rbuige — 1, indicating that the profiles reach to the re- 
gion dominated by disk light according to the bulge-disk 
decomposition. 

Three effects may influence the shape of the minor axis 
velocity dispersion profiles in our galaxies, (i) the intrinsic 
velocity structure of the bulge, (ii) the galaxy inclination, 
which determines which projection of the velocity ellipsoid 
is sampled by the data, and (iii) the brightness of the disk 
in the region sampled by the kinematic profiles. 

Because most profiles reach the region where the sur- 
face brightness is dominated by the disk (Fig.JHJ), flat dis- 
persion profiles imply similar dispersions for bulge and 
host disk. At the median inclination of the sample, i = 
73°, minor-axis dispersion profiles preferentially sample 
the radial component of the velocity dispersion ellipsoid. 
Our flat profiles then indicate a surprisingly high radial ve- 
locity dispersions in the regions dominated by disk light, 
and an equally surprising continuity of the radial velocity 
dispersions into the region of the bulge. Steep dispersion 




Fig. 8. Velocity dispersion profiles for our sample of galax- 
ies. Bulges with centrally-peaked velocity dispersions pro- 
files are drawn in the top panel, whereas bulges with 
flat velocity dispersion profiles are plotted in the bottom 
panel. In the abscissa, radii have been normalized to the 
minor axis bulge radius (Tabled Col. 12), defined as the 
radius along the minor axis where bulge and disk have 
equal surface brightness. 



profiles, in contrast, indicate distinct radial dispersions in 
the bulge and host disk, with dispersions increasing to- 
ward the center as is expected in dynamically hot, self- 
gravitating ellipsoids. 

The disk central surface brightness influences the ve- 
locity dispersion gradient. Fig. shows the distribution 
of V(<t) against the face-on disk central surface brightness 
/io- Bright disks have shallow dispersion profiles, while the 
steepest profiles are found in faint disks. We find it unlikely 
that this dependency arises from biases in the derivation 
of the kinematic profiles. Our simulations fSect. ET5|) show 
that dispersions as low as a ~ 50 km/s are well reproduced 
by the FOURFIT algorithm. Selection biases would ham- 
per the detection of disk dispersions for faint disks first, 
but dispersion gradients of all values are found for these 
disks. The lack of steep dispersion profiles in bright disks 
appears therefore real. 

Fig. [2b shows the distribution of V(cr) against e#(r e ), 
the ellipticity measured at the bulge effective radius. A 
trend is found, in the sense that more flattened bulges 
show shallower dispersion profiles. This dependence does 
not arise from viewing-angle effects: Fig. |Hh shows that 
V(a) does not depend on the galaxy inclination. Shallow- 
profile objects may be kinematically described as thick- 
ened disks. In low-ellipticity bulges, the a profiles increase 
to ward the cen t er, ag ain as expected from hot ellipsoids. 

iKormendvl (^93) discusses that some bulges share 
strong similarities with disks, on the basis of their low ve- 
locity dispersions, high V mdiX /<Jo for their ellipticity, and 
nuclear spiral structure. The disky nature of bulges is usu- 
ally discussed in relation to late-type, low-mass bulges 
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Fig. 9. Dependency of the velocity dispersion gradient on galaxy parameters, (a) V<r against face-on i^-band disk 
central surface brightness /io,D- (b) Va against bulge ellipticity e#. (c) V<j against disk ellipticity cd- Barred or 
box-peanut shaped galaxies in our sample are represented with open triangles. 




. Our dispersion profiles indicate that 



this situation affects also early-types, including several 
SO's. 



It appears that the disky nature of bulges cannot be es- 
tablished on the basis of spheroid luminosity; our velocity 
dispersion gradients do not correlate with bulge luminosity 
or with central velocity dispersion. Gradients do not cor- 
relate either with morphological type index, bulge Sersic 
index n, bulge and disk scale lengths and bulge effective 
surface brightness. We have searched for a dependency on 
the presence of bars or box-peanut shaped bulge isophotes. 
Such galaxies are depicted with triangles in Fig. While 
bar detection is difficult at the inclinations of our sample 
galaxies, it is apparent that barred and box-peanut bulges 
do not occupy preferred regions in Fig. El 



The massive bulge of NGC 5746 (Sb) is the only ob- 
ject showing a strongly-defined central velocity disper- 
sion minimum. A weak central minimum is also found in 
NGC 5443. Central velocity dispersion minima are gen- 
erally understood as tracing embedded disks, bars, or 
low-mass spheroids. Therefore, given the abundance of 
isophotal substructure in the nuclei of these galaxies, it 
is surprising that central velocity dispersion minima are 
so infrequent among our objects. As mentioned above, 
our sample includes six box- or peanut-shaped bulges. For 
these objects, the bar producing the peanut may not reach 
the center; or, alternatively, the central dispersion may 
be affected by the presence of secondary nuclear stellar 
bars/disks, with dispersion minima in bars being found 
for given viewing angles only. In NGC 5746 the local cen- 
tral minimum may be due to the main bar reaching all 
the way to the inner arcsecond, or to the presence of a 
secondary bar or nuclear disk. 



4.3. High-order Gauss-Hermite terms 

Profiles of the hs skewness term of the LOSVD are con- 
sistent with zero in most cases. In the few cases with sig- 
nificant non-zero hs (NGC 5475, 6010, 7331), hs does not 
have an opposite sign to the mean velocity that is typical 
of major axis profiles. This is reasonable as any rotation 
detected in the spectra is due to slight slit misalignment 
respect to the minor axis. 

Similarly, the kurtosis profiles h± are generally flat 
and close to zero. In 12 galaxies however, /14 is signifi- 
cantly positive above 0.1 somewhere in the profile, denot- 
ing centrally-peaked LOSVDs. 

4.4. The CaT profiles 

We have measured the CaT* profiles for our sample of 
galaxies in order to trace correlations between kinematic 
and population features along the galaxies minor axis. Our 
analysis of the profiles does not show any correlation with 
either kinematic (V. <r, ha , I14) or photometric parameters 
ffrom iBalcells et ai1l2003|) . The lack of correlation might 
be due to the very small dependence of the CaT* on age 
and met allicity for the values we expect for this type of 
galaxies (|Vazdekis et alJl2003jl. 

In a forthcoming paper ((Falcon-Barroso et al.l 12003^ 
we will present an extensive study of the near-IR Call 
triplet - (7 relation for our sample of bulges. In that pa- 
per we give measurements of the recently defined near- 
infrared Call triplet indices (CaT, CaT*), Paschen (PaT) 
and Magnesium (Mgl) indices (see Cenarro et al.ll200l4 
for our sample of bulges. We show that both the CaT* 
and CaT indices decrease with central velocity dispersion 
a with small scatter. The decrease of CaT and CaT* with 
cr contrasts with the well-known increase of another a- 
element index, Mg2, with a. 



12 



J. Falcon- Bar roso et al.: Minor axis kinematics of 19 SO-Sbc bulges 



5. Conclusions 

We have extracted minor- axis kinematic profiles from 
spectra in the Ca II near-IR region for a well-studied sam- 
ple of early - to inter m e diate- type, inclined galaxies from 
the iBalcells fc Peletierl (jl994h sample. SSP models from 
IVazdekis et alJ ([20031) at 1.5 A resolution have been used 
in addition to the traditional template stars. We have com- 
pared results obtained with the different templates and 
determined that synthetic models produce slightly bet- 
ter results than those from stellar templates, in the Call 
triplet region. We discuss the advantages/disadvantages of 
using SSP models instead of stellar templates. Our main 
conclusion is that synthetic models, when coming from 
a well constructed stellar library, help to reduce several 
drawbacks such as the template mismatch problem af- 
fecting all the traditional kinematical methods. We find a 
good agreement with high-quality kinematic profiles from 
the literature. SSP model spectra used as templates have 
clear applications for kinematic studies of galaxies at high 
redshift. 

We find minor axis rotation in almost half of our galaxy 
sample. We discuss two types of minor axis rotation: ro- 
tation in the outer parts, probably due to slit misalign- 
ments, and inner rotation, which are associated to disky 
deviations of the isophotes from pure ellipses. We also 
find some cases of kinematically offset nuclei (NGC5707 
and NGC7331) for which we do not find a clear expla- 
nation. We show that flattened bulges tend to show shal- 
low velocity dispersion gradients and similar dispersions in 
bulge and inner disk, possibly indicating a thickened disk 
structure for these bulges. Finally, we find that our veloc- 
ity dispersion gradients correlate with disk central surface 
brightness and with bulge ellipticity, but do not correlate 
with disk ellipticity. Gradients do not correlate either with 
bulge luminosity, bulge central velocity dispersion, bulge 
Sersic index n, bulge effective surface brightness, galaxy 
morphological type, and disk or bulge effective radii. 
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Appendix A: Profiles 

In this section we present kinematic profiles for each 
galaxy (heliocentric velocity, velocity dispersion, I13, I14), 
toge ther with an i7 - band NIC MPS (F160W) contour map 
from lPeletier et al.1 (|l999h (except for NGC 7332 for which 
a WFPC1 image, taken from the archive, is plotted in- 
stead). The length of the 1.2 arcsec wide slit delimits the 
range along the minor axis profile from which spectra has 
been binned to obtain a S/N>20. The CaT* index (see 
Sec. 13. 6|) is shown in the bottom panel. Filled circles rep- 
resent data points with signal to noise ratio above 33, 
whereas filled triangles represent data points with signal 
to noise between 20 and 33. Contour levels have 0.5 mag 
difference. The left hand side of the profiles, for all galax- 
ies, corresponds the dust-free side, while the right hand 
side corresponds to the dusty side. The PA (on top of 
the figure) is the position angle (N-E) of the dust-free mi- 
nor axis. Error bars determined in the simulation for each 
S/N (see sec. 13.5)1 are located on the right hand side of 
each panel, S/N=33 (short bar) and S/N=20 (long bar). 
A dotted line has been drawn, on each panel, at the level 
of the central value for the velocity, velocity dispersion and 
CaT* index and it is at the zero value for the hs and 
parameters. Vertical dotted lines indicate the minor axis 
bulge radius and center of the galaxy (see table [l] for more 
details). 
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